Further Indications of Jet Rotation in New Ultraviolet and 

Optical HST/STIS Spectra 1 



Deirdre Coffey 1 , Francesca Bacciotti 1 , Thomas P. Ray 2 , Jochen Eisloffel 3 , Jens Woitas 3 



o 
o 



(N 

m 
O 

o 

a: 
6 

CO 

c3 



ABSTRACT 

We present survey results which suggest rotation signatures at the base of T Tauri jets. Ob- 
servations were conducted with the Hubble Space Telescope Imaging Spectrograph at optical and 
near ultraviolet wavelengths (NUV). Results are presented for the approaching jet from DGTau, 
CWTau, HH30 and the bipolar jet from TH28. Systematic asymmetries in Doppler shift were 
detected across the jet, within 100 AU from the star. At optical wavelengths, radial velocity 
differences were typically 10 to 25 (±5)kms _1 , while differences in the NUV range were consis- 
tently lower at typically 10 (±5)kms _1 . Results are interpreted as possible rotation signatures. 
Importantly, there is agreement between the optical and NUV results for DGTau. Under the 
assumption of steady magnetocentrifugal acceleration, the survey results lead to estimates for the 
distance of the jet footpoint from the star, and give values consistent with earlier studies. In the 
case of DG Tau, for example, we see that the higher velocity component appears to be launched 
from a distance of 0.2 to 0.5 AU from the star along the disk plane, while the lower velocity com- 
ponent appears to trace a wider part of the jet launched from as far as 1.9 AU. The results for 
the other targets are similar. Therefore, if indeed the detected Doppler gradients trace rotation 
within the jet then, under the assumption of steady MHD ejection, the derived footpoint radii 
support the existence of magnetized disk winds. However, since we do not resolved the innermost 
layers of the flow, we cannot exclude the possibility that there also exists an X-wind or stellar 
wind component. 
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Collimated jets a s sociated with yo ung stars 
(|Eisloffel et al. 2000t iBallv et al. 20071) are be- 
lieved to play a key role in the star formation 
process, as they may be able to extract a substan- 
tial amount of the excess angular momentum from 



the star /disk sy stem (e.g. Rav fc Bacciotti 2003t 
iRav et al. 2007 ). Several models propose appar- 
ently fejasflolej™^ 
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iFerreira 1997tlshu et al. 2000llPudritz et al. 2007 
IShang et al. 2007 ) but the observational evidence 
favouring one model over another has remained 
elusive. Difficulties arise because of the small 
spatial scales involved (since jet launching occurs 
within a few AU of the star), and from the fact 
that many such protostars are heavily embedded. 
However, in recent years observations of a small 
number of jets emitting in H2 lines, and of jets 
from T Tauri stars which can be traced back to 
the source, have provided new and interesting in- 
sights. 

In particular, a new parameter for constrain- 
ing theoretical models has emerged from recent 
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high angular resolution observations: indications 
of rotation within these jets about their symmetry 
axis. This interpretation followed from the detec- 
tion of systematic radial velocity differences across 
the jet in high angular resolution spectra of the jet 
borders t aken with the slit placed parallel to the 
jet axis (iDavis et al. 200d : Bacciotti et al. 2002t 
Woitas et al. 2005j T More recently, we have un- 
dertaken a survey with the Hubble Space Tele- 
scope Imaging Spectrograph (HST/STIS) of eight 
jets, each from one of six pre-main sequence stars, 
in order to determine if signatures of rotation are 
common among these objects. For this study, 
performed at optical wavelengths, we adopted an 
instrument orientation in which the slit was per- 
pendicular to the outflow direction. Relative to 
the parallel-slit orientation, this considerably re- 
duced the exposure time needed to extract a pos- 
sible rotation signature. The results for the first 
phase of observations in this survey, which in- 
cluded the bipolar jets from RWAur and TH28, 
showed gradients in Dopple r shifts across th e 
jet of 10 to 25(±5)kms -1 (|Coffev et al. 2004 . 
A number of theoretical analyses then followed 
demonstrating the implications of interpreting 
measured Doppler gradients as jet rotation. That 
is, when combined with other jet parameters, 
measured radial velocities differences allow the 
extraction of important information such as 
the jet launching radius, under the assumption 
of st eady magnetohydrody n amic (MHD) ejec- 
tion (lAnderson et al. 20031: IPesenti et al 2 004: 



Woitas et al. 2005t iFerreira et al. 20061 ). 

While the results of the first phase of our 
survey were very encouraging, the limiting res- 
olution (spatially and spectrally) of HST/ST1S 
in the optical wavelength range restricts us to 
the peripheral regions of the jet. However, 
emission from protostellar outflows is domi- 
nated by lin es in the ultraviolet . In particu- 
lar, models dHartiean et al. 1987 ) and observa- 
tions ( Hartigan et al. 19991 ) of shocks both show 
that low excitation Mg II doublet emission at 
2796 A and 2803 A is stronger in protosteller out- 
flows than 'traditionally' observed optical lines 
such as the [O I]AA6300,6363 doublet. Therefore, 
HST/STIS observations were scheduled to exam- 
ine the same targets at near ultraviolet (NUV) 
wavelengths. Three jet targets, of the five al- 
located, were observed before the failure of the 



HST/STIS power supply on August 3rd 2004. 

We report here on the results of the sec- 
ond phase of observations for our optical survey, 
and the observations of the NUV survey. We 
present a radial velocity analysis of the optical 
dataset, which constitutes three jet targets i.e. 
the approaching jets from T Tauri stars DG Tau, 
CWTau and HH30. We also present the NUV 
dataset, which also constitutes three jet targets 
but from two T Tauri stars, i.e. the approach- 
ing jet from DG Tau and the bipolar jet from 
TH28. Our observations are described in Sec- 
tion^ while the method of analysis and results 
are detailed in Section[3l In Section|4j we discuss 
the measured parameters in the context of the 
magneto-centrifugal acceleration process, and we 
summarise our conclusions in Section[5] 

2. Observations 

Spectroscopic observations, in both the optical 
and NUV wavelength regions, were made at the 
base of jets from several T Tauri stars (listed in Ta- 
bleQJ. Jet position angles were determined from 
archival HST images. The observing procedure 
used involved centering the HST/STIS slit on the 
T Tauri star, rotating the slit to a position angle 
perpendicular to the jet axis, and then offsetting 
the slit to a position along the jet which is a frac- 
tion of an arcsecond from the source, as illustrated 
in Figure[T] 

HST/STIS optical observations were conducted 
for the approaching jets from T Tauri stars 
DG Tau, CWTau and HH30 (proposal ID 9435 
- observations of the o ther jets in the sam e pro- 
posal are described in Coffev et al. 2004h . The 
three jets were observed with slit offsets of 0'.'3, 
0'.'2 and 0'.'6 respectively. These correspond to 
a deprojected distance along each jet of 68 AU, 
37 AU and 84 AU. The larger offset for HH30 
was deemed more appropriate since the star and 
the jet base are obscured by the system's edge- 
on disk. The CCD detector was used with the 
G750M grating, centred on 6581 A, and a slit of 
aperture 52x0.1 arcsec 2 . Spectral sampling was 
0.554 A pixel -1 , corresponding to a radial veloc- 
ity sampling of ~25 kms -1 pixel -1 , and spatial 
sampling was 0. "05 pixel -1 . The effective spec- 
tral resolution, achieved using Gaussian fitting, is 
typically one fifth of the sampling, while the spa- 
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Fig. 1. — Orientation of the jet and slit for each target. The arrow on the slit indicates the positive direction 
of the y-axis on the position-velocity contour plot, Figure^ The jet orientation can be compared with 
compass orientation given in the upper corners of each radial velocity profile, Figures|4] and [5] The arrow 
around the jet axis indicates the sense of rotation as implied by the results in this paper. The requested slit 
position angle was 90° with respect to the value of PA jet , TableHJ The actual slit position angles for DG Tau 
and HH 30 differ from the requested values by +3° and -6° respectively. This was due to problems during 
observations in finding the right combination of guide stars for the requested angle. The sense of rotation is 
denoted according to the viewpoint of the observer looking down the approaching jet towards the star. 
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Target 


Location 


Distance 


M* 


V S ys 




PA jet 


References 






(pc) 


(Mq) 


(kms -1 ) 


(deg) 


(dog) 




DGTau 


Taurus 


140 


0.67 


+ 16.5 


52 


226 


1) 2 


CW Tau 


Taurus 


140 


1.4 


+ 14.5 


41 


155 


3, 4, 5 


HH30 


Taurus 


140 


0.45 


+21.5 


1 


33 


6, 7, 8 


TH 28 


Lupus 3 


170 




+5 


10 


98 


9, 10 



Table 1: Details of TTauri jet targets investigated in this paper, 
after correction for the systemic heliocentric radial velocity, v s 



All radial velocity results, v ra d, are quoted 
ys . The inclination angle of the jet, z jet , 



is given with respect to the plane of the sky. Values for the jet position angle, PA jet , were determined 
from archival HST images. For our observations, we req uested a slit position angl e of 90° with respect to 
PA jek Figured! Reference s - ( I) lEisloffel et al. 19981 C2 l iBacciotti et al. 20021: C3)lG6mez de Castro 19931 
f4l iHartmannet al 1981 (5) [Hartigan et al 2004 ; (6) Petv et al. 20 061: (7) lAp penzeller et al. 20051: (8) 
Mundt et al. 19901: (9)lGraham fc Hever 19881: (ToT Krautter 1986 . 



tial resolution is twice the sampling (i.e. 2 pixel 
resolution of O'.'l). One long exposure of ^2100 s 
was made of the approaching jet from DG Tau on 
Dec 1 2003 and the approaching jet from CWTau 
on Dec 21 2003, while two long exposures (to be 
co-added) were made of the approaching jet from 
HH30 on Nov 30 2003. The data were calibrated 
through the standard HST pipeline, subtraction 
of any reflected stellar continuum was performed, 
and signals from defective pixels were removed. 

HST/STIS NUV observations were conducted 
for the approaching jet from DG Tau, and each 
lobe of the bipolar jet from TH28 (proposal ID 
9807). A slit offset of 0."3 was used in each 
case. For TH28, this corresponds to a depro- 
jected distance along the jet of 52 AU. The NUV 
MAMA detector was used with the E230M echelle 
grating, centred on 2707 A, and a long slit of 
aperture 6x0.2arcsec 2 to ensure the full width 
of the jet was observed. Spectral sampling was 
0.045 A pixel -1 , corresponding to a radial velocity 
sampling of ^5 kms -1 pixel -1 , and spatial sam- 
pling was 0. "029 pixel -1 . The effective spectral 
resolution achieved was on average 5 kms -1 , while 
the spatial resolution was twice the sampling (i.e. 
2 pixel resolution of 0'.'058). Two long exposures 
of ~2500s (to be co-added) were made of the ap- 
proaching jet of DGTau on Dec 1 2003, and each 
lobe of the bipolar jet from TH 28 on Jun 18 2004. 
Although the data were largely processed through 
the standard HST pipeline, the combination of an 
echelle grating with a long slit meant that wave- 
length calibration was not conducted as part of the 
pipeline procedure. Wavelength calibration was 
therefore necessary in the initial stages of data re- 
duction, and was carried out using standard IRAF 



routines. 

3. Results 

In our chosen observing mode, the slit does 
not cover the star. The instrument spatial 
line spread function has a half-width at zero- 
maximum of 0'.'15 in the optical regime and 0'.'145 
in the NUV. Therefore, observed emission lines 
in both the optical and NUV regions originate in 
j et material excited t hrough low-velocity shocks 
(jHartigan et al. 19871 ). All optical targets were 
found to emit in the forbidden emission line 
doublets [O I], [N II] and [S II] (i.e. at vac- 
uum rest wavelengths of 6302.046, 6365.537 A, 
6549.86, 6585.27 A and 6718.29, 6732.67 A respec- 
tively). All three jet targets observed in the NUV 
were found to be strong emitters in both lines of 
the Mg II doublet (i.e. at vacuum rest wavelengths 
of 2796.352 A and 2803.531 A). Their profiles are 
broad, extending over some 200 kms -1 , consistent 
with the idea of shock excitation. For each tar- 
get, the extension of the broad line profile in the 
NUV is approximately the same as for the optical 
lines presented in this and previous studies but, 
spatially, the emission region appears to be more 
concentrated on the jet axis. In the NUV, spa- 
tial extention is represented by FWHM values in 
the range of 0'.'12 to 0'.'20 for both jet lobes from 
TH 28 (as opposed to 0'.'27 and 0'.'41 for [O I] and 
[S II] receding jet lines respectively - from optical 
spectra taken with HST/STIS, proposal ID 9435). 
Values for the approaching jet from DG Tau are 
similar, with FWHM values of O'.'l 2 in the NUV 
(as opposed to 0'.'15 and 0'.'20 for [O I] and [S II] 
optical lines respectively). The trend shows that 
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the high excitation NUV lines preferentially orig- 
inate closer to the axial region of the jet. This is 
consistent with the expectations of jet models with 
shock heating, according to which layers closer to 
the axis are more excited than those further away, 
since the velocity of the shock front is higher on- 
axis and gives rise to higher shock temperatures 
there. 

3.1. Qualitative analysis 

If a jet is rotating, we would expect to ob- 
serve a difference in the magnitude of the Doppler 
shift between material on one side and the other 



of a jet's symmetry axis (see ICoffev et al. 2004) . 
Therefore, signatures of jet rotation are expected 
to present themselves in the form of a contour tilt 
in position-velocity diagrams of jet spectral emis- 
sion lines. In other words, the tilt in contours illus- 
trates that radial velocities are lower on one side 
of the jet, thus giving a qualitative indication of 
systematic radial velocity differences as expected 
from a rotating jet. 

Figures[2] and [3] show position- velocity contour 
diagrams for selected optical and NUV emission 
lines, the contour levels for which are given in Ta- 
ble[2] In each contour plot, the positive direction 
of the y-axis corresponds to the slit direction as 
illustrated in Figured) While a position-velocity 
contour tilt can generally be recognised at optical 
wavelengths, a tilt is not so apparent in the NUV 
range. It is likely that the broad profile shape in 
the dispersion direction combined with the rela- 
tively narrow spatial FWHM may detract from a 
visible delineation of any trend. The existence of 
an absorption feature for the NUV permitted tran- 
sitions may also make any subtle contour trend less 
evident. Nevertheless, quantitative measurements 
yielded positive results. 

3.2. Quantitative analysis 

Quantitative measurements of jet rotation re- 
quire determination of differences in radial veloc- 
ity at positions equidistant either side of the jet 
axis. The intensity peak in the spatial direction of 
the jet emission is assumed to mark the physical 
jet axis. The physical jet axis must be aligned with 
the centre of a pixel in order to measure equidis- 
tant either side of it. In practice, emission peaks 
were offset slightly (in most cases by about 0.25 



pixels from the central pixel row) due to offsets in 
telescope and instrument pointing. Therefore, the 
first step in the analysis was recentre the emission. 
The method was to perform a Gaussian fit on each 
column (i.c in the spatial direction), identify the 
column with the highest amplitude fit, determine 
the offset of its peak, and recentre the 2D emission 
according to this offset amount using an interpola- 
tion technique. The precision reached was ± O'.'Ol 
(i.e. one fifth of the spatial sampling of 0'.'05). 
Offsets were double-checked by binning the emis- 
sion onto the spatial axis and Gaussian fitting the 
peak. The offsets were in agreement for the opti- 
cal data, but only the latter technique was relied 
upon when processing the NUV data. Note that 
the position-velocity diagrams (Figures[2] and [3]) 
are shown prior to this recentering, specifically to 
illustrate the presence of contour tilts before the 
data were processed in any way. 

Gaussian fitting was then used to determine 
where the emission peaked in each pixel row par- 
allel to the jet axis (i.e. in the dispersion direc- 
tion). From this analysis, radial velocities were 
determined at each distance from the axis and a 
radial velocity profile of the jet in the transverse 
direction was obtained. Radial velocity profiles are 
plotted in Figures|4] and [5] (in which the orienta- 
tion, derived from Figure[TJ is given in the upper 
corners of each plot). Where the low and high 
velocity jet material is clearly distinguishable, the 
plots have been divided in to low and high velocity 
components. Plots reveal that the on-axis jet ma- 
terial is travelling at the highest velocity, while the 
borders of the jet are travelling at lower velocities. 
However, the velocity profiles are not symmetrical. 
Velocities on one side of the axis are lower than on 
the other. Furthermore, each emission line traces 
a different degree of asymmetry. In general, [O I] 
lines are clear tracers of radial velocity asymme- 
try while [S II] lines show much more symmetric 
profiles. This combination of contributions from 
different emission lines and different velocity com- 
ponents can draw from the clarity of some of these 
profile plots. The clearest case is the plot for the 
TH28 receding jet at optical wavelengths. 

The differences in radial velocities were cross- 
checked by comparing them to the results of a 
cross-correlation routine. While the latter has the 
advantage of being independent of the line pro- 
file shape, it was found to be very sensitive to 
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Fig. 2. — Position-velocity contour diagrams in selected emission lines for each jet targets observed in the 
optical regime. Note that the [O I] emission traces both the higher and lower velocity jet material, while 
the [N II] emission is only seen at high velocities and the [S II] emission is mainly concentrated at low 
velocities. Furthermore, higher velocity emission appears to be more concentrated on the jet axis, while 
lower velocity emission is spatially extended. The skew in contours is indicative of jet rotation. The positive 
direction of the y-axis is illustrated pictorially as the slit direction in Figure[T] Plots are corrected for 
systemic heliocentric radial velocity, v sys , TablefTJ Contour values of the fluxes for each panel, in units of 
erg cm -2 sec -1 A -1 arcsec -2 , are given in Table[2j 



Target 


Emission 


Contour 


Contour 


Contour 




line 


Floor 


Ccilin 




Interval 


DG Tau approaching jet 


[O I]]A6300 


1.0xl0" ib 


4.4x10" 


13 


2.0X10" 14 




[N II]A6583 


5.0xl0" 15 


8.1x10" 


14 


4.0X10" 15 




[S II]A6716 


3.0xl0" 15 


7.1x10" 


13 


4.0X10" 15 




Mgll A2796 


l.OxlO" 13 


2.1x10" 


12 


l.OxlO" 13 


CW Tau approaching jet 


[O I]A6300 


3.0xl0" 15 


3.2x10" 


14 


l.OxlO" 15 




[N II]A6583 


2.0xl0" 15 


1.3x10" 


14 


7.0X10" 16 


HH 30 approaching jet 


[O I]A6300 


7.0xl0" 15 


6.4x10" 


14 


3.0X10" 15 


TH 28 approaching jet 


Mgll A2796 


l.OxlO" 13 


3.9x10" 


13 


2.1X10" 14 


TH28 receding jet 


Mgll A2796 


2.0xl0" 13 


2.0x10" 


12 


1.3xl0~ 13 



Table 2: Contour intervals of the emission line surface brightness of position-velocity diagrams, Figure^ in 
units of erg cm" 2 s" 1 A" 1 arcsec -2 . 
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Fig. 3. — Position- velocity contour diagrams in Mg II A2796 emission for the jet targets observed in the NUV 
regime. The positive direction of the y-axis is illustrated pictorially as the slit direction in Figure[T] Plots 
are corrected for systemic heliocentric radial velocity, v sys , Table[TJ Contour values of the fluxes for each 
panel, in units of ergcm~ 2 sec -1 A^ 1 arcsec -2 , are given in Tabled 
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the extent and level of the background. Although 
there is generally good agreement between the 
methods (within 1 or 2kms~ 1 ), Gaussian fitting 
proved more reliable in the optical region where 
all emission lines were of suitable profile shape. 
In the NUV range, in the case of DG Tau where 
the higher velocity material escaped absorption, 
Gaussian fitting proved suitable and was found to 
be in reasonable agreement with cross-correlation 
(within 2kms _1 ). In the remaining cases, Gaus- 
sian fitting was not reliable due to absorption 
of part of the profile. In these cases, measure- 
ment of radial velocity differences relied on cross- 
correlation. Tables[3] and 0] give radial velocity dif- 
ferences, Av ra d, for individual emission lines in 
each jet lobe. These measurements are also illus- 
trated graphically in Figure[6] Asymmetries were 
typically 10 to 25kms _1 in the optical region, and 
lOkms- 1 in the NUV region. 

Figure[7] is shown as an illustration of the ef- 
fect of a decrease in signal-to-noise on the three- 
sigma error bars associated with each Gaussian 
fit using IDL. Data points below signal-to-noise of 
3 were generally deemed unreliable. Given suffi- 
cient signal-to-noise, we are able to routinely mea- 
sure radial velocity differences to one fifth of the 
velocity sampling, conservatively speaking. This 
yields an accuracy of ±5kms _1 in the optical and 
±lkms _1 in the NUV. However, since additional 
uncertainties are introduced by the NUV absorp- 
tion feature, we prefer to give the conservative er- 
ror estimate of ±5kms _1 for our NUV results. 

Meaningful analysis of much fainter NUV emis- 
sion lines evident in the observations was not pos- 
sible. One of these emission peaks was identi- 
fied, in the DG Tau jet spectrum, as the combined 
emission lines of C II] AA2324,2325,2326,2327,2328 
travelling at 180 km s -1 . The other emission peak 
appears to be [O II] A2470 , again if a jet radial 
velocity of about 180 km s^ 1 is assumed. 

The Mg II absorption feature dominates the jet 
emission from TH28, more so than in the case 
of DG Tau. This becomes evident when emission 
is binned for plotting in one dimension, Figure[8] 
The ratio of Mg II fluxes integrated across each jet 
for the entire emission line is given in each plot. 
Fluxes are quoted without extinction correction, 
which could not be determined from these spectra 
because pairs of transitions, from the same up- 
per level, for any given species are of insufficient 



separation in wavelength to allow extinction cal- 
culation. Also, no compensation has been made 
for the absorption feature since the exact emission 
profile shape is not known. The absorption dip is 
located at low blue-shifted velocities for all three 
targets. Approximate speeds were measured as - 
26kms~ 1 in the TH28 receding jet, -15kms _1 in 
the TH28 approaching jet and -43kms _1 in the 
DG Tau approaching jet (i.e. with repect to the 
vacuum wavelength but with no adjustment for 
the systemic velocity). 

4. Discussion 

4.1. Indications of Jet Rotation 

In almost all cases, the results show clear gra- 
dients in Doppler shift across the jet base, in the 
optical and NUV wavelength regions. One inter- 
pretion of these radial velocity differences is as in- 
dications of jet rotation close to the source where 
the jet is launched. 

Position-velocity diagrams for the optical for- 
bidden emission lines give a qualitative indication 
of rotation in the form of a tilt in the contours, Fig- 
ureO Furthermore, from these diagrams and the 
NUV flux profile plots (Figure[5]) we see how differ- 
ent emission lines trace different velocities within 
the jet: [N II] lines appear to trace only higher ve- 
locities; Mg II and [O I] lines trace both higher and 
lower velocities; and [S II] lines trace mainly lower 
velocities. The diagrams also show how the jet 
becomes spatially resolved at low velocities, as is 
clearly evident in cases where the emission traces 
both higher and lower velocity material, i.e. the 
DG Tau [O I] and [S II] emission and the CWTau 
[O I] emission. Meanwhile, for all targets, the [S II] 
lines show this spatial extension but little or no 
rotation, in the sense that radial velocity differ- 
ences mainly lie within the error bars of 5kms _1 
about zero (Figure[6|). This may be due to the 
fact that at a given distance from the jet axis, 
line-of-sight and projection effects (as discussed in 
(Pesentiet al. 2004)) reduce the measured Doppler 
shift more considerably in the case of the broader 
[S II] flow than, for example, the [O I] flow. 

The transverse radial velocity profiles for all 
targets (Figures^] and O further highlight the 
asymmetry in the radial velocities across the jet. 
The clearest cases are those of the DG Tau ap- 
proaching jet and the TH28 receding jet. Here 
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Fig. 4. — Radial velocity profiles across the jets in various optical emission lines. All radial velocities are 
corrected for systemic heliocentric radial velocity. Profile asymmetries were interpreted as indications of jet 
rotation. 
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Fig. 5. — Radial velocity profiles across the jets in NUV emission lines (similar to Figure|4|). The optical 
datasets are overlaid on the NUV datasets, to illustrate how wel l the results for the two wavelength regions 
fit together. Optical data for TH 28 were published previously (| Coffey et al. 2004). The NUV profiles for 
the TH 28 bipolar jet should be taken as an indication only, because Gaussian fitting proved very difficult 
and so an accurate profile could not be constructed. The radial velocity differences in Tabled were obtained 
by relying on a cross-correlation routine, as described in Section]! 
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Target Distance A.v ra d for An rQC j for Av ro£ j for Ai; rac f for A.v ra d for At? rac f for Av ro£ j for 

from axis [O I]A6300 [O I]A6363 [N II]A6548 [N II]A6583 He IA6678 [S II]A6716 [S II]A6731 





(arcsec) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


DGTau 
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14 


7 


2 


5 


4 


approaching 
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10 (-2) 
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17 (-3) 
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CWTau 
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15 (11) 
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1 (3) 
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15 (0) 






15 






-5 (4) 


jet 
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12 (5) 












... (-1) 


HH30 


0.05 
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-1 


-1 


approaching 


0.10 
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3 




-1 


-2 


jet 


0.15 


4 










-3 


-3 




0.20 


1 












-4 




0.25 














-2 



Table 3: Radial velocity differences, Av rac i, between one side of the jet axis and the other measured using 
single Gaussian fitting. Only where the emission is divided into a clear higher and lower velocity component 
was double Gaussian fitting used. In these cases, the results for the lower velocity component are given 
in brackets beside the results for the higher velocity component. Where blanks appear in the table, the 
emission was either Doppler shifted off the CCD, or was too faint to decipher. The Ha line from the jet 
was contaminated by reflected emission from the star thus rendering it too difficult to analyse. Finally, 
the [O I]A6363 line was often contaminated by the signal from a defective pixel. The error bars on the 
data analysis are at most ±5 km s _1 but reduce with increasing signal-to-noise (see Figure[7]). A signal-to- 
noise cut-off ratio of 3 was used, as measurements below this level were generally deemed unreliable, which 
therefore excludes the other regions of the flow. 



Target Distance Av ra d for AiVad for 

from jet axis Mg II A2796 Mg II A2803 





(arcsec) 


(kms -1 ) 


(km s 


DGTau 


0.029 


7 


2 


approaching 


0.058 


14 


6 


jet 


0.087 


14 


11 




0.116 


8 




TH28 


0.029 





-3 


approaching 


0.058 


2 




jet 


0.087 


1 




TH28 


0.029 


-3 


"1 


receding 


0.058 


-3 


-2 


jet 


0.087 





-4 




0.116 


1 


-3 




0.145 


1 


-2 




0.174 


3 






0.203 


6 






0.232 


3 





Table 4: Radial velocity differences, Av ra d, between one side of the jet axis and the other, measured in NUV 
lines using a cross-correlation technique. Cross-correlation was found to be more suitable than Gaussian 
fitting given that the profiles include disruptive absorption. Since the velocity profiles of Figure[5] rely on 
Gaussian fitting, they should not be expected to exactly match this table, but are intended only as an 
indication of the transverse velocity profile. Where blanks appear in the table, the emission was too faint 
to decipher. The accuracy reached with the data analysis was approximately ±5 km s . A signal-to-noisc 
cut-off ratio of 3 was used, as measurements below this level were generally deemed unreliable (Figure[7]), 
which therefore excludes the other regions of the flow. 
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Fig. 6. — Transverse radial velocity differences measured in optical and NUV lines, for a slit placed across 
the jet less than 100 AU from the star. These values are taken from Tables[3] and [4] with the exception of 
those for the optical lines in the TH 28 receding jet, which are taken from Coffey et al. (2004). The horizontal 
strip marks the maximum error bar in the velocity difference about zero, but is more representive of the 
typical error on NUV data points than optical data points (see Figure[7]). 
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Fig. 7. Left - Velocity profiles for selected target and emission lines plotted with error bars representing 
the three-sigma error associated with each Gaussian fit; Centre - signal-to-noisc of binned intensity profile 
after emission is centered on the jet axis; Right - variation of signal-to-noise with errors on Gaussian fits. 
Fitting becomes unreliable typically below a signal-to-noise ratio of 3. 
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Fig. 8. — Flux profiles of the Mg II AA2796,2803 emission doublet, integrated across the jet, for the approach- 
ing jet from DG Tau and each lobe of the bipolar jet from TH 28. The velocity scale is corrected for systemic 
heliocentric radial velocity, v sys ^ Table[TJ Fluxes are quoted without compensation for the absorption feature, 
as the exact emission profile shape is not known. 
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the plots reveal that the on-axis jet material is 
travelling fastest while the borders of the jet are 
travelling slower. The result is a v-shaped pro- 
file. But these profiles are not symmetrical: ve- 
locities on one side of the axis are lower than on 
the other, a natural outcome if the jet is rotat- 
ing. In the case of the CWTau higher velocity jet 
profile, the emission is more confined. Its profile 
does not incorporate a significant border contri- 
bution, and it is clearly separated from lower ve- 
locity material (Figure[2]). It therefore does not 
trace an asymmetric v-shape, but instead exhibits 
a simple slope whereby velocities on one side of 0" 
are lower than on the other side. It is natural in 
such cases that the central datapoint should not be 
the highest radial velocity point, but should have 
an intermediate velocity value (within error bars). 
Examining all target profile plots, and allowing for 
an uncertaintly of 5kms _1 , the variety within the 
plots can be explained according to one or both 
of these two basic profile shapes. The combined 
profiles of the optical and NUV datasets illustrate 
how measured radial velocities are comparable for 
Mg II and optical lines for both the approaching 
jet from DG Tau and the receding jet from TH 28. 
The approaching jet from TH 28 does not follow 
this trend. The fact that this jet is inclined at 
only 10° with respect to the plane of the sky (Ta- 
bleQJ means that it has low radial velocity val- 
ues, and so the peak jet velocity almost coincides 
with the velocity of the absorption feature. This 
prevents reliable estimates of radial velocity from 
NUV lines for this target. 

Measured radial velocity differences, Tables[3] 
and 21 show values typically in the range of 
10 to 25kms- 1 for DGTau and CWTau in 
the optical region within 0."25 of the jet axis. 
Meanwhile, marginal differences were detected for 
HH30 of typically Skms^ 1 , which lie close to or 
within the error bars about zero. In the NUV 
region, the Mg II doublet shows radial veloc- 
ity differences for DGTau of about 5 to 15 km s -1 
within 0."1 of the jet axis which is in line with, 
but notably lower than, the optical values. The 
values for TH28 are all within error bars about 
zero, except for a measurement of 6kms _1 at 0."2 
from the jet axis in the receding jet. Note that 
in a previously published paper, optical values of 
10 to 24kms~ 1 wer e measured for the same region 
(jCoffev et al. 20041 ). Difficulties in tackling the 



absorption dip possibly contribute to lack of de- 
tectable velocity differences in this NUV dataset. 

To comment generally on Tables[3] and [4j there 
appears to be a drop in velocity differences close 
to the jet axis for DG Tau in the NUV. Such lower 
values close to the axis were also noted for other 
jets in the optical analysis of Coffey et al. (2004) 
and outlined as being most l ikely due to line-of- 
sight and projection effects (jPesenti et al. 2004; 
Dougados et al. 20041 ). However, this drop-off 



does not appear in the optical results for DG Tau 
to the same extent. Also, Coffey et al. (2004) 
describes [N II] lines as showing higher rotation 
since this emission line traces higher velocities, but 
this is not supported here in the optical results of 
DGTau and CWTau, or the NUV results where 
the Mg II doublet traces equally high velocities as 
[N II] emission. 

Interpreting these gradients in Doppler shift as 
indications of jet rotation reveals a sense of rota- 
tion as indicated in Figure[T] Supporting this in- 
terpretation is the fact that there is agreement in 
direction of Doppler gradient between the DG Tau 
optical and NUV results. For the TH 28 receding 
jet, there is no clear gradient detected in the NUV, 
although the outlier of 6 kms -1 at 0'.'2 would sug- 
gest the same sense of rotation as measured in 
previously publishe d optical data for this target 
(|Coffev et al. 20041 ). 

A large velocity difference was not measured for 
HH30. This is surprising, given that this system 
is in the plane of the sky and so should present the 
most easily measureable transverse Doppler gradi- 
ent. Interpreting the differences as rotation would 
imply that the approaching jet is rotating clock- 
wise, looking down the approaching jet towards 
the star. This finding would conflict with th e 
measured sense of disk rotation (|Petv et al. 2006t) . 
However, since these measurements are marginal, 
it appears that we cannot conclusively derive a 
sense of rotation for the jet. Meanwhile, no ro- 
tation signature was detectable in the surround- 
i ng CO outflow associated with this optical jet 
(|Petv et al. 2006h . 

An alternative explanation for a gradient in ra- 
dial velocity across the jet base could be the pres- 
ence of asymmetric shocks, since shocks are in- 
voked to explain jet line emission. Indeed it is 
conceivable that we are observing a combination 
of effects, including both jet rotation and asym- 
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metric shocking of jet material. However, obser- 
vations of Doppler gradients across the bipolar jet 
of RW A ur taken with sever al slits parallel to the 
jet axis (jWoitas et al. 2005) demonstrate that if 
an asymmetric shock was causing these gradients, 
the shock velocity asymmetry would have to per- 
sist along the jet for at least ~ 200 AU to explain 
these observations. Not only that, but the same 
shock asymmetry would have to persist over this 
distance in both lobes simultaneously. Further ob- 
servations of Dopp ler gradients in bipolar systems 
(jCoffev et al. 20041 ) . which also reveal agreement 
in the direction of the gradient within the bipolar 
jet, demonsatrate that the direction and magni- 
tude of the shock velocity asymmetry would have 
to be the same in each lobe in the bipolar system. 
Therefore, while asymmetric shocking is a conceiv- 
able explanation, it becomes less likely within the 
context of these results. 

Finally, Cerqueira et al. (2006) have shown how 
the small-angle precession of a jet can produce the 
transverse velocity differences seen here and hence 
mimic rotation. In particular, they have modelled 
a jet using a mean velocity of 300 km s -1 along 
with a precessional angle and period of 5° and 
8 years respectively. Some of the jet parameters 
are based on the DGTau jet but, for example, 
they choose a jet inclination angle of 45° rather 
than the observed v alue for DGTau of 37.7°±2 
(|Eisloffel et al. 1998h . For their simulations they 
also assume that the jet has a time-dependent ejec- 
tion velocity that varies with the same period as 
the precession. In the case of DG Tau, there is 
ample evidence from its morphology that the jet 
is precessing (Lavalley-Fouquet et al. 2000; Douga- 
dos et al. 2000). Thus for DGTau a precessional 
model may be reasonable; although this does not, 
of course, exclude the presence of rotation as well. 
That such a model is applicable however in every 
case in which we see transverse velocity differences 
of a few tens of kms -1 seems contrived. For ex- 
ample, if the precessional angle were only 1-2°, 
instead of 5° , no measurable effect would be seen 
in their pure precessional model of the DG Tau jet. 
Likewise the precessional period, etc., are critical. 
Moreover, as Cerqueira et al. (2006) remark them- 
selves, there is no evidence in some cases, e.g. the 
RWAur jet, for precession. That said, any am- 
biguity as to the cause of the transverse veloc- 
ity differences in the case of individual outflows 



can, and should, be removed through modelling. 
In particular precessional parameters can be ob- 
tained from multi-epoch imaging and fitting of jet 
emission centroids as a function of distance from 
the source. 

4.2. Error Analysis 

With reference to errors, it should be mentioned 
that, during the first stage of data reduction (see 
the beginning of Section [3~2|) , when each emission 
peak was recentred to the middle of the pixel in 
the central row of the detector (i.e. the nominal 
0"position on the detector), a trend was observed 
in the pixel offsets from the 0"position for each 
emission line. Peak offsets ranged from 0.1 to 1.5 
pixels depending on the target, although offsets in 
emission lines for a given target varied by <0.5 
pixels. This highlighted three types of possible in- 
strument misalignments, which had the potential 
to contribute to position-velocity contour tilt and 
thus mimic rotation. The first is physical tilt (i.e. 
a tilt of the slit with respect to the nominal ob- 
serving position angle), and has values of +0.22 
degrees in the optical region and +0.27 degrees 
in the NUV. The second is optical distortion tilt 
(i.e. the slit image on the detector is curved and 
tilted) which varies in angle across the detector 
depending on the grating used, and has values of 
+0.001 degrees in the optical region and +0.9 to 
+1.2 degrees in the NUV. The third is position an- 
gle error (i.e. a difference between the slit position 
angle used in the observations with respect to that 
requested), and could be on the order of degrees. 
For example, as noted in Figured] the actual slit 
position angles for DG Tau and HH 30 differ from 
the requested values by +3° and -6° respectively. 
This was due to problems during observations in 
finding the right combination of guide stars for the 
requested slit position angle. 

All possible corrections for errors arising from 
instrument misalignment were made before data 
analysis was conducted. For the optical data, the 
HST pipeline calibration accounted for physical 
tilt. The optical distortion, on the other hand, was 
not removed in the pipeline. Moving from right to 
left along the CCD, the spectrum gradually curves 
upwards with a total difference of half a pixel from 
one end to the other. In other words, the [O I] 
lines are raised higher on the detector than the 
[S II] lines. This was adequately accounted for 
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by recentering the emission peaks in the middle of 
a pixel. There could still remains elements of the 
distortion which cannot be accounted for, since its 
exact nature is unknown, but the effect on mea- 
surements cannot be significant. A confirmation 
of this is given by HST/STIS archival data with 
almost the same instrument configuration as our 
observations. A long-slit spectrum of Mz 3 was ob- 
tained with HST/STIS on June 23, 2002 (Proposal 
ID 9050). The same grating (G750M) and detec- 
tor (CCD) were used, but the slit aperture size 
used was slightly different (i.e. 52x0.05arcsec 2 
rather than 52x0.1 arcsec 2 ). These slits are in fact 
coaligned and their optical paths are the same, 
so this difference in aperture size does not affect 
our investigation. The archival spectrum shows 
flux levels comparable in intensity and spatial ex- 
tension to our data, and exhibits [S III] emission 
which hits the detector in the region between the 
two components of the [O I] doublet. Reassur- 
ingly, no tilt is present in the [S III] line. Further- 
more, Gaussian fitting to pixel rows either side 
of the 0"position for the [S III] line revealed veloc- 
ity differences of effectively zero (ie . < . 5 km s ~ 1 ) . 
Given the extreme unlikelihood that instrumental 
effects could operate in this manner (i.e. causing 
a tilt in the [O I] doublet but no tilt in emission 
between the doublet lines), this data confirms the 
fact that the instrument is not introducing spuri- 
ous contour tilts in our data. 

For the NUV data, standard IRAF routines 
were used to wavelength calibrate the data. These 
routines rectify the spectrum both spatial and 
spectrally, although without accounting for the op- 
tical distortion. Therefore, prior to this rectifica- 
tion, the calibration lamp lines were used to deter- 
mine the extent of the optical distortion tilt in the 
region on the detector of interest. The STIS cal- 
ibration reference tables were then manually ad- 
justed to incorporate the value of this tilt angle 
prior to wavelength calibration. In this way, stan- 
dard routines accounted for both physical tilt and 
optical distortion. As in the case of the optical 
data, any residual effect cannot be accounted for 
since the nature of the distortion is unknown. Nev- 
ertheless, it cannot significantly affect our mea- 
surements since, for DG Tau, our NUV measure- 
ments are in reasonable agreement with the optical 
measurements even though they were taken with 
different instruments and have a different optical 



distortion to contend with. 

The remaining tilt to be addressed is possible 
inaccuracy in the position angle specified for ob- 
servations. Although this is mainly corrected for 
through recentering the emission peaks to the cen- 
tre of a pixel (as described in the beginning of 
Section l3~2|) . it is possible that a tilt error could 
remain unaccounted for. The precise effect of this 
misalignment is not known, since the kinematic 
nature of the underlying flow is not known. How- 
ever, in the simplest case of a cylindrical, non- 
accelerating, non-rotating flow a slit misalignment 
does not introduce a radial velocity gradient across 
the jet. Therefore, for our observations, we con- 
sider the effect insignificant. 

The fact should be noted that we see the same 
sense of Doppler gradient with respect to the slit 
orientation for all T Tauri jets in the optical and 
NUV survey for which a Doppler gradient is de- 
tected. This includes previously published opti- 
cal results for both lo bes of the bipolar jets from 
TH28 and RWAur (ICoffev et al. 20041 ). (Note 
that the previously published results for LkHa 321 
have since proved inconclusive.) This fact might 
be regarded as unusual, and could point to instru- 
mental error. Analysis of archival Mz3 [S III] emis- 
sion (as described above) excludes this possibility. 
Furthermore, bare in mind that, of the eight jet 
lobes, the seven that show Doppler gradients (i.e. 
excluding LkHa321) are each launched from one 
of five T Tauri systems. This effectively reduces 
the sample size, making the statistic more accept- 
able. It should be noted that a rotation analysis of 
the radial velocities within the jets from T Tauri 
stars DG Tau (Bacciotti et al. 2002) and RWAur 
(Woitas et al. 2004) conducted in the optical re- 
gion with the slit placed parallel to the jet axis 
yielded results in agreement with our observing 
mode in which the slit is perpendicular to the jet 
axis. This strongly negates instrumental error. To 
further strengthen the case, it is important to also 
highlight the fact that the observations for the jets 
from both DGTau and TH28 in the optical and 
NUV regions were conducted with different detec- 
tors (and so different optics). For DGTau the re- 
sults are in good agreement, while for TH28 the 
results are not inconsistent. This combination of 
studies, along with the careful analysis of instru- 
mental effects, strongly suggests that, although it 
seems fortuitous that all jets (with the possible 
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exception of LkHa321) have Doppler gradients in 
the same direction with respect to the slit orienta- 
tion, it is not impossible and does indeed seem 
to be the case. Finally, note that four of the 
five T Tauri systems with Doppler gradients are 
located in the same cloud (i.e. Taurus- Auriga) , al- 
though no particular pattern of system alignment 
could be readily identified. 

4.3. Mg II Absorption 

In recent studies, several spectra of T Tauri 
stars have been obtained, in which the as- 
sociated jet was included in the slit coverage 
(|G6mez de Castro fc Verdugo 200ll ; lArdila et al. 20021) . 
Ardilaet al. (2002), in particular, present a spec- 
trum of DGTau taken with the Goddard High 
Resolution Spectrograph onboard HST. The cov- 
erage was ~2" centered on the star, and so the 
spectrum included both the magnetospheric ac- 
cretion region and the base of the jet. The au- 
thors assume that all the Mg II emission originates 
from the magnetosphere, and attribute the broad 
blue-shifted dip to absorption by a cool outflow. 
Figure[5] shows a comparison of their data and 
ours, and clearly demonstrates that the jet is both 
absorbing and emitting. Their spectrum includes 
emission of Mg II radiation from the jet such that 
our spectrum corresponds both in velocity and 
intensity to the blue-shifted plateau reported in 
Ardilaet al. (2002). In fact, one-dimensional NUV 
spectra similar to ours were obtained for the bow- 
shock HH47 A with the Faint Obj ect Camera on- 
board HST (jHartigan et al. 19991 ). In that case, 
the line intensities have been successfully modelled 
as arising from composite shocks advancing into 
a low density, and moderately ionised, medium. 
Figure[9] also serves to illustrate that the stellar 
PSF does not enter the slit in our observations. 
Recall from Section[3] that, for our observations 
of the jet base, the wings of the HST spatial line 
spread function in both the NUV and the optical 
regions do not extend as far as the star, and so 
there can be no contamination from photospheric 
line emission. This is supported in Figured] by 
the lack of emission at intermediate velocities in 
our observations compared to observations cen- 
tered on the star. The higher jet velocities (> 
250 kms -1 ) of Ardilaet al. (2002) are not seen in 
our spectrum, which seems to indicate that our 
observations were taken at a distance from the 



star which is before the end of the jet acceleration 
zone. This is consistent with evidence of increas- 
ing velocity with distance along the DG Tau jet as 
reported by Bacciottiet al. (2000). Similarily, in- 
termediate velocity emission (~100kms _1 ) either 
originates from the wings of the stellar emission 
profile, or from jet material closer to the star which 
is yet to be fully accelerated. 

In our spectra, a low velocity absorption feature 
is also present within the emission profile of all 
three targets (Figure[8]). This dip is travelling at 
low blue-shifted velocities, has a narrower FWHM 
(~40kms _1 ), and is not significantly asymmetric 
as would be expected from wind absorption. Since 
it is 6/we-shifted in all the observed targets, includ- 
ing the receding jet from TH 28, it cannot originate 
in the jet lobe itself. The dip does coincides with 
interstellar absorption velocities, and so it is con- 
ceivable that it originates from interstellar cloud 
material. Mg II interstellar absorption features 
are present, for example, in spectra of main se- 
quence F stars (Bohm-Vitense et al. 2001). If this 
were the case in our spectra, explanation of the ab- 
sorption FWHM would necessitate the combined 
contribution of two or three interstellar clouds 
(see also discussion by Ardilaet al. (2002)). This 
assumption of multiple interstellar clouds would 
have to be made in the direction of both DG Tau 
and TH28, making this interpretation less plaus- 
able. The feature could be more readily explained 
by an expanding shell around the T Tauri system. 
In contrast to attributing it to a cooler outflow 
wind, an expanding shell would be observed as a 
blue-shifted dip regardless of whether the observed 
jet is approaching or receding. A shell interpre- 
tation would also account for a FWHM which is 
possibly too large to be caused by ISM absorption. 

A final point concerns the flux intensity ra- 
tio Mg IIA2796/A2803 within the jet. Statisti- 
cally, this ratio is expected to be 2:1 when emit- 
ted by an optically thin gas, while it approaches 
1:1 in conditions of opt ical thickness (see e.g. 
Bohm-Vitense et al. 2001 ). Hence, the flux ra- 



tios (Figured]) for the TH 28 receding jet and the 
DG Tau approaching jet demonstrate the gas to be 
optically thin. Meanwhile, the TH 28 approaching 
jet is evidently optically thick, and the FWHM 
of the dip is slightly larger here than for the other 
targets. Therefore, in this case, the absorption fea- 
ture is very likely to incorporate a self-absorption 
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component. Besides, this jet is known to be heav- 
ily embedded in reflection nebulosity (Graham & 
Heyer 1998). This combined contribution to the 
dip, along with the fact that the emission and ab- 
sorption peaks overlap considerably, leads to par- 
ticular difficulties in determining accurate radial 
velocities for emission from jet target. 

4.4. Jet Launch Point 

Under the assumption that the radial velocity 
gradient measurements described in Section[3] are 
to be interpreted soley as jet rotation, and that 
the mechanism for launching this portion of the jet 
can be accurately described by a magnetocentrifu- 
gal stationary wind, then the quantities derived in 
this work can be used to constrain the properties 
of the launch model. In particular, we can use our 
radial velocity measurements, in combination with 
the values in Table[TJ to And a range for the radius 
from the star on the disk plane of the jet foot- 
point (or launch point), zuq. This was calculated 
using Equation[T] in Anderson etal. (2003), which 
is derived from general conservation properties of 
magnetocentrifugal models for jet acceleration, 

x ( v ^ r 4/3 fiM 1/3 a) 

VlOOfcms- 1 / \IM Q J y ' 

where vJoo is the radius from the jet axis of 
the observation (conducted at an effectively infi- 
nite distance from the star), t>0 iOO is the toroidal 
velocity at the observation radius, v p .oa is the 
poloidal velocity at the observation radius, and 
M* is the mass of the star. Since the inclination 
angle, ij et , given in Table[T]is with respect to the 
plane of the sky, then = (Aw rac ;/2)/ cos i and 
v p = (vrad)/ sin i, where v ra d is the average of the 
radial velocity measured at a given distance each 
side of the jet axis. The equation is valid for val- 
ues of Vcj, 3> v p , which is achieved at distances from 
the star of typically a few AU, and hence well be- 
fore the position of our slit at several tens of AU. 
Where M* is unknown, it is assumed to be equal 
to 1M Q . 

The range for the jet footpoint radius for each 
target is given in TableO Values were calculated 
for a sample of emission lines, rather than giv- 
ing an average for each target, since each line 



traces different levels of radial velocity. Calcula- 
tions were carried out for radial velocity differ- 
ences measured far from the jet axis, to minimise 
line-of-sight and projection effects. These effects 
only allow measurement of lower limits in radial 
velocity within typically one spatial FWHM from 
the axis, a critical point for the determination of 
true toroidal velocities and hence footpoints (Pe- 
sentietal. 2004). 

From the results of the three emission lines for 
DG Tau, for example, we see that the higher ve- 
locity component appears to be launched from a 
distance of 0.2 to 0.5 AU from the star along the 
disk plane, while the lower velocity [O I] com- 
ponent appears to trace a wider part of the jet 
launched from as far as 1.9 AU. These values are in 
the same range as those estimated for the DG Tau 
jet of 1.8 AU by Bacciottiet al. (2002), using veloc- 
ity differences from data taken with the STIS slit 
parallel to the jet axis, and of 0.3 to 4 AU by An- 
derson et al. (2003) for the same dataset (where the 
analysis was carried out on only the lower velocity 
component of the jet). The results for the other 
targets also fall into this range, determining the 
footpoint of the jet to be within 0.2 to 3.9 AU from 
the star. Note, also, that DG Tau presents a good 
case for examining the variation of poloidal veloc- 
ity with footpoint radius. Here the low and high 
velocity jet material have comparable radial veloc- 
ity differences, and so we can see that decreasing 
the poloidal velocities gives rise to an equivalent 
increase in footpoint radius. 

Given the uncertainties associated with many 
of the physical parameters (such as stellar mass, 
inclination angle and stellar distance), these cal- 
culations represent estimates of the jet footpoint 
and so no error bars are given. Nevertheless, the 
results support the idea that jets are launched, via 
the magneto-centrifugal mechanism, at footpoint 
radii within a few AU of the star. The results con- 
sistently show that, given our spatial resolution, 
the portion of the jet that we observe appears to 
be launched from a region up to a few AU from 
the star on the disk plane. However, this one-to- 
one mapping of the footpoint is probably an over 
simplification. In fact, MHD instabilities are likely 
to influence the geometry of the field lines beyond 
the Alfven surface, and so the jet material can- 
not necessarily be so easily traced back to exact 
footpoints on the disk plane. 
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Target 


Emission 
line 


(arcscc) 


(AU) 


(km s- 1 ) 


Vrad 

(km s _1 ) 


(km s~ x ) 


(km s _1 ) 


T*70 

(AU) 


DG Tau approaching jet 


A6300 


0.15 


21 


26 (18) 


195 (60) 


21 (15) 


248 (76) 


0.5 (1.9) 




A6583 


0.15 


21 


7 


165 


6 


209 


0.3 




A2796 


0.116 


16 


8 


193 


6 


244 


0.2 


CWTau approaching jet 


A6300 


0.15 


21 


12 


108 


8 


164 


0.6 




A6583 


0.10 


14 


15 


102 


10 


155 


0.5 


TH 28 receding jet 


A6300 


0.20 


34 


16 


8 


8 


46 


3.9 




A6583 


0.20 


34 


24 


27 


12 


155 


1.0 




A2796 


0.20 


34 


6 


13 


3 


75 


1.0 



Table 5: The radius from the star on the disk plane of the jet footpoint (or launch point), wq, calculated 
for targets in the optical and NUV using the method described in Anderson et al. (2003). Values in brackets 
relate to the lower velocity component. The mean radial velocity, taken from values equidistant either side of 
the jet axis, v ra d, is quoted as an absolute value (and after systemic heliocentric radial velocity correction). 
Toroidal and poloidal velocities were calculated as v$ = (Av ra d/2)/ cos i a nd v p = (v ra d)/ sin i. The TH28 
[O I] and [N II] results are taken from data reported in (jCoffev et al. 20041 ). 
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Fig. 9.— Mg IIA2796 flux profile for the DG 
Tau system. The solid line represents our obser- 
vation of the approaching jet at 0'.'3 from the star, 
taken from Figure[8] The dashed line represents 
observations in 1996 of a 1.74 arcsec squared area 
centered on the star, extracted from Figure 5 of 
The comparison serves to il- 



Ardila et al. 2002 
lustrate that the outflow is both absorbing and 
emitting in Mg II, and also that the stellar PSF 
does not enter the slit in our observations. 



In the context of a HST/STIS survey to ascer- 
tain the commonality of T Tauri jet rotation, we 
report on observations of three jet targets at op- 
tical wavelengths (i.e. the approaching jet from 
DG Tau, CWTau and HH30), and three jet tar- 
gets at NUV wavelengths (i.e. the bipolar jet from 
TH 28 and the approaching jet from DG Tau). Al- 
most all targets show distinct and systematic ra- 
dial velocity asymmetries for opposing positions 
with respect to the jet axis, within 100 AU from 
the source. Radial velocity differences of typically 
10 to 25 (±5) km s _1 were found at optical wave- 
lengths, while values at NUV wavelengths were 
lower at typically 10(±5)kms~ 1 . We interpret 
these radial velocity asymmetries as rotation sig- 
natures in the region close to the star where the 
jet has been collimated but has not yet manifestly 
interacted with the environment. Possible instru- 
ment contribution to error has been thoroughly 
examined and ruled out. We also note the pres- 
ence of a Mg II blue-shifted absorption feature 
in all three NUV cases. This possibly originates 
from the combined contribution of several inter- 
stellar clouds, or from an expanding shell of ma- 
terial around the T Tauri system. 

It was hoped that the higher spatial and spec- 
tral resolution of the HST/STIS in the NUV would 
afford a more quantitative analysis, in terms of an 
error bar reduction from one fifth of the spectral 
sampling in the optical region (i.e. ±5kms _1 ) 
to that of the NUV (i.e. ilkrns" 1 ). However, 
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we did not expect to find the tricky combination 
of a broad profile shape in the dispersion direc- 
tion, a relatively narrow spatial FWHM, and sig- 
nificant interruption by absorption. This forced 
us to adopt a conservative approach to errors, i.e 
±5kms -1 . Nevertheless, in the case of DGTau, 
the NUV measurements fit well with the results 
obtained from the optical data, and there is agree- 
ment in magnitude and direction of the Doppler 
gradients within the two wavelength ranges. The 
results for TH28 were, unfortunately, less clear 
and no sense of Doppler gradient could be estab- 
lished with certainty. Finally, the survey results 
are consi stent with previously p ublished results for 
DGTau (|Bacciotti et al. 20021 ). 

The determined rotational velocities were 
shown to be in agreement with magneto-centrifugal 
launch model predictions (e.g. Anderson et al. 2003) 
In the case of DG Tau, for example, we see that the 
higher velocity component appears to be launched 
from a distance of 0.2 to 0.5 AU from the star 
along the disk plane. Meanwhile the lower veloc- 
ity component appears to trace a wider part of the 
jet launched from as far as 1.9 AU. The results for 
the other targets are similar and consistently show 
that, at the resolution of the observations, the jet 
is launched from a region of up to 3.9 AU from the 
star on the disk plane. Therefore, if indeed the 
detected Doppler gradients trace rotation within 
the jet then, under the assumption of steady MHD 
ejection, the derived footpoint radii support the 
existence of magnetized disk winds. However, 
since we do not resolved the innermost layers of 
the flow, we cannot exclud e the possibility t hat 
there also exists an X-wind dShang et al. 20071) or 
stellar wind component (jPudritz et al. 200 7). Al- 
though we cannot probe the inner axial part of the 
flow at current resolution, this may be possible in 
the near future via interferometry. 
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